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Calculations on Countercurrent Electromigration ! 

By G. Breit 2 and F. L. Friedman 



The mathematical problems presented by the countercurrent electrolysis method for 
isotope separation are treated under the simplified assumptions described in the text as 
(a), (b), (c), (d), (e), and (f) in section II, 2, of this paper. The solutions are worked out 
in section II, 3, for a tube of finite length. The formulas are put into numerical form in 
equations 9e, 9f, etc. Simpler formulas are obtained in section II, 4, in the approximation 
of tubes of infinite length. Some of the results are presented in the form of graphs that are 
described in section III. The orders of magnitude of effects to be expected are discussed in 
section IV. 



I. Introduction 

The countercurrent electromigration metho< I 
has been found to give isotope separation in the 
case of potassium [1, 2] 3 and chlorine [3]. In 
the present report some diffusion problems of this 
method are treated mathematically. Formulas 
are given for the calculation of the concentration 
along the tube. Some typical numerical examples 
are given to illustrate the use of the formulas and 
graphs and to show the variation of concentration 
with time and distance. 

II. Calculations 

1. Notation 

7ii=number of more mobile ions per cubic 

centimeter (supposedly light). 
n h = number of less mobile ions per cubic 
centimeter (supposedly heavy). 
lii. fjL h = mobilities of light and heavy ions respec- 
tively. 
D h D h = effective diffusion coefficients of light 
and heavy ions. 
v= velocity of liquid flow (positive when 
toward cathode). 



1 This paper will appear in volume 6, division III, of the Manhattan Project 
Technical Series. 

2 Professor of Physics, University of Wisconsin, now, Yale University. 
Prof. Breit Was on the staff: of the National Bureau of Standards in connec- 
i ion with the initiation of the atomic-energy project when this work was done. 

3 Figures in brackets indicate the literature references at the end of this 
paper. 



isoelectric field (positive when directed 
from anode to cathode). 
n h nl= normal values of n h n h . 

M = common symbol for m or fi h if distin- 
guishing between them is immaterial. 
Z)=common symbol for T) x or B h if distin- 
guishing between them is immaterial. 
x= distance along tube measured from 

anode to cathode. 
V= volume of cathode compartment. 
A = cross soetional area of tube free to carry 

electric current. 
An = fx l — \x h . 
L= length of tube. 

t=time counted from beginning of experi- 
ment. 
L' = V/A= effective length of cathode. 

2. Simplifying Assumptions and Basic Equations 

It is assumed that: (a) The mobilities and the 
diffusion coefficients have the same values through- 
out the tube. This is not true, but it is felt that 
the variation of these quantities may be neglected 
for the present preliminary estimates, (b) The 
electric field is assumed to be constant along the 
tube. This is also not strictly true. As the 
electric field could be maintained at a constant 
value, this assumption should not matter much. 
(c) The cathode is assumed to be stirred per- 
fectly, (d) The concentrations in the anode are 
supposed to be maintained at their normal 
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values, (e) The solution is supposed to be com- 
pletely dissociated, (f ) Forces acting on ions due 
to space charge effects are neglected. 

At a point in the tube having distance, x, from 
the anode the flux of light ions is 






(1) 



The first term takes account of the combined 
action of the electric field which gives rise to the 
migration and of the counter current velocity, v. 
One has, therefore, 



bt 



bn t 

~bi 



+| [n,G.,*+.) -D,3g]=P (2) 

f + l!i \n*<J»E+v)-D h ^f]=0 (2a) 



for the equations expressing equilibrium of the 
light and heavy ions within the tube. At the 
cathode compartment x=L. Since the liquid 
inside the cathode is stirred one has 

^^=^ l =^ l E+v)n l -D l ^(x=L) (2b) 



Vdn, 



=<£„= (» h E+v)n h -D h %g (x=L) (2c 



It is conyenient to introduce 



a t =niE+v, a h =n h E+v 




(2d) 


a] at 




(2e) 






(2f) 


\ — a% T \ — a " T 




(2g) 


V a, V . L' 
AD,~AL* l -L 1 "'' y "~ 


L' 

L X " 


(2h) 



This amounts to introducing units that are some- 
what more natural for the discussion than the 
original units. In terms of the natural units one 
has instead of eq 2 and 2a; 



bni__ b 2 n t bn t bn h _ b 2 n h bn h 
*rr b& &£ (*tj~ b& bi h 



(3) 



The conditions for the cathode compartment 
become 

/ bnA ( bn\ ( bn h \ ( bn h \ 

(3a) 

The problem is defined by equations 3, 3a, and the 
conditions 

ni(o,t)=n h n h (o,t)=n° h (3b) 



n l (x,o)=n° h n h (x,o)=n h 



(3c) 



which express the fact that at all times the con- 
centrations are normal at the anode and also that 
at t=0 the concentrations are normal throughout 
the tube. 

Since the equations are of the same form for nt 
and n h , the indices I and h may be dropped in 
most of the work. 

3. Tube of Finite Length 

In this section an exact solution for a tube of 
finite length is given. The motion of one isotope 
only is considered and the formulas are supposed 
to be applied by working out separately the con- 
centrations of the two isotopes. The indices I, h 
are omitted in the discussion of the behaviour of a 
single isotope. The equation 



bn_ bn , b 2 n 
b~r~~ ~^^W 



(4) 



is being solved subject to boundary conditions (eq 
3a, 3b, and 3c). 

For r=co one has in addition 



n(£, co)=n°et. 



(4a) 



The solution can be constructed by expanding the 
initial state in terms of a complete set of functions. 
For 7=0 the functions are mutually orthogonal 
and the method is then well known. For y^O the 
functions are not mutually orthogonal. It is 
nevertheless possible to give an explicit usable 
solution in this case also. 

Equation 4 is first transformed by 



n=eV*N 



It assumes the form 



dW N__bN 

d£ 2 4 br 



(4b) 



(4c) 
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The function n is expanded as 

»(£ T )=n eS+e! /2 2 aMi)*"** ( 4d ) 

where t\(£)e~ 5,r are solutions of eq 4c so that 

3gr+(/J.- j)*.=0. (4e) 

The functions t?< are taken so as to satisfy the 
boundary conditions 

»<(0) -o, [^] { , x =(^+|) y '( x ) ( 4f ) 

The first of these corresponds to eq 3b, the second 
to eq 3a. The condition in eq 3c is taken care of 
by making 



are introduced. In terms of these, eq 5c gives 



2n°siiih!+S<W<tt)=0 



(4g) 



All initial conditions are satisfied automatically 
as a result of eq 4f and 4g. It remains to deter- 
mine the a u 

Combining two previous equations one has 

p i (X)^(X)-©;(X)»i(X)=(ft-j8 l ) Ptv(0f>,(0<ft 

Eliminating v' if v\ by means of eq 4f, one has for i^j 

J»i(0»itt)«+7»«(X)^(X)=O. (5) 

In order to expand the function 

-2n° sinh (f/2) =/ft)=Sa^({) (5a) 

the functions # t -(£) are first normalized to unity: 

J\(£)dS=l (5b) 

By integration one obtains from eq 5a, 

at-yJ>tMv,(\)a,=ft 9 ft= PviWWt (5c) 
The abbreviations 

0,=0«(X), cr=Sa^=/(X) (5d) 



a* 



' 1+7^ 



(5e) 



The boundary conditions in eq 4f and the nor- 
malization in eq 5b are satisfied by 



,„. V2/Xsmh(rf/2) ,- — 

^o(€)=-M= ^ w \ t ; r= V 1 -40, r=real number 
ysinh rX/rX— 1 

(6) 
(6a) 



where 



tanh ^r=^ 



2 l+ T (l-r 2 )/2 
determines r. 

Equation 4c is similar in form to the time de- 
pendent wave equation in wave mechanics. Only 
real 0^0 need be considered, therefore. Ac- 
cordingly, only real or pure imaginary values of 
r are of interest. 

For 7^0, X>0 eq 6a has in addition to r=o at 
most one real root for r and this lies in 0</<l. 
In addition there are pure imaginary roots 



r=ip s 



These satisfy 

fcan (^") = r 



-7(l + p/)/2 



(6b) 



(6c) 



There is an infinite number of roots p n forming a 
discrete set. 

If the velocity of the liquid is higher than the 
migration velocity then X and y have negative 
values. From eq 6a one obtains 



1 r i+rl + 7 (l-r)/2 "| 
> - r m Ll-r l + 7(l+r)/2j 



(6d) 



For 7< —2, values of r in 0<><1 + 2/y give 

4/| 7 



values of X in 



l-2/|7| 



+ 7 (l+r)/2. 
in 0<r<] 
>\>-oo. For t> -2 



there are no real roots for r. For a<^0 there is 
thus also at most one real root for r. 

Performing the calculation according to eq 4d, 
5a, 5e, and making use of eq 6a and 6c, it is found 
that 



8 f*v* 



sinh -y sinh -^- 6 



x- 1+ ^+^ 



2 V 2)y 1 + 2) b*-l) 

0+1)-*? 



(7) 
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where the summation covers real and imaginary sinh (r s X/2) is replaced by i sin (p s X/2). 
r s . For imaginary r s eq. 6b and 6c are used and As a special case 



-id -Or 



l(\,r)=nV-tW°S79 
A s = o( & 



2 ' 



gives the concentration at the cathode. For large 
positive 7 the real root moves toward unity so 
that r =l — e where e<<0- I n this limit one has 
approximately 



G-'^.)[(>+i)'-M-K>+l)<'-^ 



"6 X -1 



(8) 



The imaginary roots contribute terms in I/7 2 and 
become small for 7-^00. Equation 7a gives in 
this limit 



n(\, T)=n°e x —n°(e* — l)e 



7(^-1) 



(8a) 



so that the approximate time constant in units 
D/a 2 is 



T° = y(e*-1), 



(8b) 



which corresponds to a time constant in seconds 
V 



If in addition X<<1 then 

VL_VL 
1 = AD~ D 



(8c) 



(8d) 



If a<0, eq 8 to 8d still hold for large values of— 7. 

For a small cathode compartment, one has the 
simpler formulas: 
n(«,r) = 

ot 8n° ^^ ^h T smh^ g (g) 

7ie " X ' £T -l+r?+2/X 



o .2 -(l""rJ)r/4 

X ^(l-r?)(- 



•1+rf+J 



-2/X) 
(9a) 
For 7=0, X— 2, the real root is in the process of 
becoming imaginary, and the corresponding term 
in the above formulas needs special consideration 
becoming indeterminate. Making X=2 + e with 
e^O it is found that: 

~ I ^-*~ , ^ sin p 4 sin (p,{/2)<r**| 



71 e 



1 



(p s /2) 2 



(9b) 



co A 

-" 1 _1_ «2 * 

1 1-rP. 



-/3.rl 



(7a) 



(9c) 



In this case the roots are known and are tabulated 
in Jahncke-Emde [4]. Substitution of numbers 
gives 

nQi, T)ln o = e 2 -6e-e° r -0.l8876e-^- 

0.06592e-fc'— 0.03336e"^ T — 0.02011e-^ T — 

0.013446-^-0.00962^-^-0.007226"^- 

0.00562e-^-0.0044966-^-0.00368e"^- 

0.00307e-"» r -0.002595e-^ T -0.00223e-^ T - 

0.0019296-^ 4r -0.001688e-^-0.0014896-^- 

0.001324«-^r-o.0011846-^ 8r -0.001066e-^ T - 
0.000965€-^-0.000877e^ 2 ^-0.00080l6-^- 

. . . (9d) 
Also 

7?° 

w(l, r)=n°6— ^{3e-^-0.15087e"^- 

-y/e 

0.04388e-^ T +0.02473e-^ r + 

0.013746-^-0.00979e"^ r - 
0.006646"^+ 0.00522e"^+ 
0.00390^-^-0.003236"^- . . . } (9e) 



The /3's in eq 9d and in 9e are given by: 

— * -- * l ft- 



1 



0.18876 



0.06592 



0.03336 



In order to see how the cathode volume affects 
the results, the above solution can he compared 
with that for X=2, 7=15.2195. The choice of y 
in this calculation is made so as to give r =0.98. 
Approximate values of the p's are px = 3.18, p 2 = 
6.30, p 3 =9.44. For the cathode, one obtains by 
means of eq 7a, 






= 7.3891 - 6.38856- - 0099 ' - 
0.000436- 2 - 78 ^-0.000046- 10 - 2t - . . . 

which shows that in this case only the first term 
having the time constant T=10lD/a 2 is of real 
importance. The approximation in eq 8b gives 
T=97D/a 2 and is suitable for an approximate esti- 
mate for this 7 and X. Because of the more rapid 
convergence of the series, calculations with large 
7 are easier than with small 7. 
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Another example of the application of eq 7a is 
\=2, r =0.90, 7 = 2.6995 for which Pl =3.331, 
p 2 = 6.394, and 

(-o) =7.3891 -6.3766- 0047 ^-0.011326- 3025 ^- 

0.001 15<r 10 - 47 '-0.00025<T 22 - 9T - . . . (9f) 
An example for negative a may also be given: 

X=~-l,-r =0.50,7=~8.111,pi = 6.359,P2=12.605, 

which gives 

(^) = - 36788 + °- 6319 ^+ 

0.00028<r 10 - 36 '+ . . . (9g) 

A consideration of these examples shows that 
for a rather large range of values of y, the term in 

— — \ 

e 4 is the last important one. The relation 
t° = y(e x — 1) for its time constant can be improved 
on by using a graphical or numerical solution of 



tanh 



rX 



and setting 



2 l + 7(l-r 2 )/2 
4 



' i— r5 

4. Tubes of Infinite Length 



(9h) 



Substitution of different r into eq 9d, 9e shows 
that the concentration at the cathode builds up 
much more rapidly than in the middle of the tube. 
For times appreciably smaller than the time con- 
stant the effect of the anode is, therefore, not 
important, and the calculation of exact solutions 
involves a needless amount of work. The initial 
stages of the building up of the concentration 
can be studied by supposing that the tube is 
infinitely long because the solution for an infinite 
tube gives practically zero concentration at the 
anode. 

In discussing a tube of infinite length it is con- 
venient to measure lengths from the cathode 
introducing 

i?=X-£ (10) 



The differential equation is 



d 2 n dn_dn 
drj 2 dr/ - dr 



(10a) 
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and the boundary condition at the cathode is 
conveniently expressed in terms of 



for which 






(10b) 



dn 



(»S-L-(*£-L- o*> 



The quantity m is the flux density of ions divided 
by a. It has proved convenient to use m in con- 

structing the solution because m—y >— has the 

simple initial condition of being 1 throughout the 
length of the tube. It is thus possible to construct 
first the solution for 7 = 0, to equate the m ob- 
tained for 7 = torn— 7 j- and derive m for any 

7. From m, one obtains n by quadratures. 
For a^>0, one finds 



l — 71° -y/r/n 



(^+r)2 



"1 + 7 



■*[-($)] 



-[- 



?^[-<^)> 



(2+7)7 

2(l + 7) 2 



V+r, t 



where 



$(x) 



V^ Jo 



+ T'-*(^+t)] <n) 



e-^dr, $(-x) = -$(x). (11a) 



For ?/=0 i. e. for the cathode the above formula 
is somewhat shorter: 



( n-n° \ 



T+l+7" 



1 + 7 



•I' I 



2(1 + 7) ~\2,) + 



2 



-t/4 



+ r 



7(7+2) 



(Hb) 



1+ T ' '"2(1 + 7) 2(l+ 7 ) 2 
and for 7 = 0, the following simple form obtains: 
n-n° ,-r -<s+» If, a A+A"|4. 



n 



<£2 

2 



(T -,+ l)[l-*(|^)] (11 



c) 
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The values of $(x) are available in Jahnke-Emde 
[4] in good detail and numerical calculation by 
means of eq. 11, lib, and lie is straightforward. 
A convenient form of m is 



m 



1 . 1 f(^+r)/2Vr 

= o+-f= e-*dx- 

* V 71 " J° 

p-rj roo 

^=7— e~*'dx+ 



■"? e (v+r)/y+T/y 2 



V^(i+t) 



x 'dx. 



r)/2Vr+^- (Hd) 



For a<0, X, 77, 7 are negative and r is positive. 
In place of eq 11, lib, lie, and lid one has 



n° 1 + t 2 



?[S+i^]['+K^)] 

»^|>+*(i^f)] 



+ 
+ 



(12) 



+ 






and for 7=0 



(12a) 



/- _k±i2_ 2 



c: (r _, +1) [ 1+# (|=i)] (12b) 



while for arbitrary 7<0 

2 Vxjo 



^=^-4= f (,+T)/2V V^- 



g-*7 r(v-r)/2-y/r 

-i=z r e~ x2 dx+ 

Vtt(1 + 7)J — 

-74fr^^ 7 7 ^ 2 ^ ( 12c ) 

V7r(l + 7J J-co 



The verification of the above equations by 
substitution into the differential equation and the 
boundary conditions is straightforward and is not 
reproduced in this report. 

III. Description of Figures 

In figure 1 there is shown the dependence of 
(n—ri) Jn° on r at the cathode for a tube of infinite 
length. The curve for zero cathode volume 
(7=0) is parabolic in shape at the origin. The 
initial rate at which the concentration builds 
itself up is seen to be very high. This is because 
the cathode has zero volume so that the initial 
rate is infinite. As time goes on the material con- 
centrated at the end of the tube diffuses through a 
volume proportional to -y^. The rate at which the 
concentrated material is being fed to this volume 
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Figure 1. 

X=oo, n=0; x are points from X=2, ->=0, 4=2. 



1. Parabola 0.17Vr'0.0200; 2. 7=0; 3. osculating parnbola 2-y/rJir; 4. tangent; 
5 7=0.1; 6. 7 = 1 

is proportional to r and the concentration is 
therefore proportional to t/Vt=V^- This is, of 
course, not an exact consideration, because the 
rate at which the concentrated material is being 
fed to the end of the tube depends on the concen- 
tration at neighboring places. The parabolic 
approximation is not bad, however, as is seen in 
figure 1, where the dotted curve above the curve 
for 7=0 is a parabola 0.17VV/.0200, chosen so as 
to give a good overall fit and the dotted curve 
below the curve for 7=0 is the osculating parabola 
2-yjT/w. The latter curve makes a second-degree 
contact with the curve for 7=0. In figure 1 there 
are shown also, as crosses, points that correspond 
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to \=2, £ = 2, 7 = 0. Those are seen to fall on the 
curve \=<x>, % = 2, 7=0 to within the accuracy of 
the drawing. A close agreement of the two curves 
can be expected because the diffusion distance in 
natural units (D/a) is of the order Vr^O.15 for 
the larger r in figure 1, while the length of the 
tube in the same units is 2. The time in figure 1 
may thus be considered so short that the concen- 
trated solution has not yet diffused to the anode. 
It is satisfactory in this case to make computations 
of the initial stage of the enrichment on the sim- 
plified assumption of an infinitely long tube 

(X=oo). 

Figure 2 corresponds to negative a i. e. to the 
removal of ions from the cathode end. In this 
case n—n° is negative and, therefore, —(n—n )/n° 
is plotted as ordinate. A large ordinate corre- 
sponds to a small concentration in this figure and 
the cathode volume is zero. For r=l, the con- 
centration is reduced to about 1/3 of its original 
value. The lower curve is drawn on an expanded 
time scale so as to show more detail in the initial 
value. 



























1 
















































**• 


^- 


-— -""* 














/ / 






















/ 

/ 

I 























3 4 5 

.3 4 .5 



X = 



Figure 2. 
oo, y=0, 71=0. ]. Remaining concentration divided by n°=nln° 



In figure 3 a comparison is made 4 between the 
relative enrichment obtained by having the 
countercurrent balance exactly one or the other 
ion type that it is desired to separate. In this and 
some succeeding figures a positive a that corre- 
sponds to feeding an isotope into the cathode is 
referred to by the subscript "pos" and a negative 
a that corresponds to washing an isotope out of 
the cathode is similarly referred to by the subscript 
"neg". The low T er curve in figure 3 corresponds 
to having the countercurrent of proper strength 
to balance the migration velocity of the slower 
isotope. If the initial concentrations of the more 
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Figure 3. 



X= oo, 7=0, 71=0. 



1 



( n— n° \ 



-1: 2. 



1 + 



/ n-n° \ 



mobile and less mobile isotopes are n°, n° h respec- 
tively, then at the time t this ratio is 



n K := n i\' 1 1 / n—n °\ "1 



(13) 



and the factor by which this concentration is 
increased is unity, plus the quantity plotted for 
the lower curve. 

If, on the other hand, the countercurrent is 
adjusted to balance the more mobile isotope, then 
the isotope of smaller migration velocity is being 
drawn out from the cathode and after a time / 



Ui 

n,~ 



n t 






-i+ 



\ ^ / neg I J 



(13a) 



The quantity in square brackets is plotted as the 
upper curve in figure 3. It represents, as before, 
the ratio of concentrations divided by its initial 
value and minus 1 : 



nijrin 

o / o~ 

njn h 



■1. 



It will be noted that for quite a while there is no 
marked difference in the two ways of adjusting the 
counter flow that have been considered. 

In figure 4 five curves are shown. All of them 
are for a tube of infinite length. The lower pair 
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of curves is for |7| = 0.1. In this pair the upper 
curve corresponds to an adjustment of counter 
flow corresponding to eq 13a, the lower to an ad- 
justment corresponding to eq 13. The less mobile 
isotope is being pumped out of the cathode for the 
upper curve, the more mobile isotope is being 
pumped into the cathode for the lower. 




/ n-n° \ 
' \ n° /po 8 '~ 



"A »° /ne, 



The three upper curves in figure 4 are for a tube 
of infinite length and zero cathode volume. 
Among these the highest curve, H, corresponds 
to countercurrent adjustment in eq 13a, the 
lowest to eq 13. The curve in between, curve M, 
corresponds to a countercurrent adjustment in 
which the arithmetic mean of the migration 
velocities of the more mobile and less mobile 
isotopes is equal to the velocity of the liquid. 

The value of a for eq 13a is — |jiij— iJL h \E &nd for 
eq 13 it is + |ju*— n h \E. In the computation for the 
middle curve in the upper set of three in figure 4 
the value of a is ±\ni—n h \E/2 and in accordance 
with eq 2e, the value of r for a given t is ){ of what it 
is for a= db \ni—Hh\E. For this reason the abscissae 
for curve M were made to be 4r. In this way the 
same value of the abscissae for the curves H, M, 
L corresponds to the same time, t. 

The curve, M, was computed by plotting 
[(^M°)pos/(W 7l °)neg]— 1 as ordinate and 4r as 
abscissae. It is seen that the curve M falls nearly 
half way between H and L. This indicates that 
the exact adjustment of the countercurrent does 
not critically affect the ratio of the abundances of 
the two isotopes as long as there is an approxi- 
mate balance of migration velocity and counter 



flow. With the idealized conditions assumed here, 
the relative enrichment in the cathode is changed 
more rapidly by having the counter flow too high 
rather than too low. 

The way in which the three typical adjustments 
of countercurrent give about the same result for 
small enrichments can be illustrated also in terms 
of a formula for small enrichments for 7=0. One 
has approximately 



W<U 8 ^l+2Vr/7r, (n/n°) neB ^l-2Vr/T 



so that 



/ n-n° \ 1 

V n° / pos ~ , ( n-n° \ 



^2Vt/tt (H,L) 



which shows the approximate equivalence of 
countercurrent adjustments, H and L. Also, for 
the countercurrent adjustments of the curve M 



l + 2Vr/7r 
l-2Vr/7r" 



li 



(n/n ), 



' (n/n°) ne 



-1^2V4t/tt, (M) 



so that the quantity M is the same function of 
4t as H, L are of r. 

For y^Oa change in the countercurrent gives 
a change in 7. For small t the three types of coun- 
tercurrent adjustment give nearly the same result. 
In this case 



(Wpos = 1 + |tc| 



'pos |To| 

where y =V(ni—n h )E/AD. Similarly 






To I 

For the arithmetic mean countercurrent adjust- 
ment, 

2r 



(n/n°) 



(n/n°) 



pos 



■li 



1+f 



|7o| 



neg 



1- 



2t_ 
|Yo| 



1 = 



4r 

ItoI 



and the result is nearly the same as if the counter- 
current were adjusted to keep one of the isotopes 
stationary with respect to the tube. 

In figure 5 the behaviour of a tube of infinite 
length is compared with that of a tube of finite 
length in the special case of zero cathode volume. 
The concentration of an isotope at the cathode 
end is plotted against r in both cases. The 
equilibrium concentration for the tube of finite 
length is taken to be e 2 =7A. Up to t=2, there is 
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Figure 5. 

7 = 0, n = 0. 1. X=oo. 2. X = 2. 

not much difference between the two tubes. The 
concentration increases to about four times the 
initial value at this r. For t>4, the slope of the 
curve for the tube of finite length decreases rapidly 
so that it will not pay to run above this value. 
The comparison of X=2 with an infinitely long 
tube has already been made for small r in figure 1. 
The initial stage of the divergence between the 
two curves can be seen in figure 6. 
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Figure 6. 

y =2, X=2, $=2. j are for X= oo. 

In figure 7 the dependence of concentration of a 
single isotope on r is illustrated for 7=0 at differ- 
ent points along a tube of infinite length. The 
migration velocity is supposed to be greater than 
the countercurrent corresponding to a positive a. 
At the cathode end the concentration increases 
steeply, being proportional to -y/r for small r. For 
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Figure 7. 

7=0. x are from X=oo; o are from X = 2. 1. w=0.00; 2. n=0.10, 3. n=0.25; 
4. n=0 50. 



points farther away from the cathode the initial 
portion of the curve is flat and is followed by an 
approximately linear rise. The initial flat stage 
can be thought of as corresponding to the time 
during which the concentration that first develops 
at the cathode diffuses to the point under con- 
sideration. In figure 7 there are also shown points 
that correspond to a tube of finite length with 
X = 2. These are seen to fall very closely on the 
curve for X= <». This is to be expected, because 
for the small r covered in this figure there is prac- 
tically no effect at the anode end so that it does 
not matter whether the concentration at the anode 
is maintained at its initial value or not. 

In figures 8 and 9 there is illustrated the varia- 



^2 " 

3 ^^**» 
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Figure 8. 
7=0. 1. r=0.16; 2. r«0.010;3, r=0.0016. 
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7 =0. I. t=0 010: 2. t=0.0016. 
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Figure 10. 



7 = l.\=oo. l.r=0.16.2. r-0.01. 



tion with distance of [(n — n )/ri ] pos i. e. of the 
fractional increase of concentration for the case 
of feeding the isotope towards the cathode. The 



cathode volume was taken to be zero for these 
figures, and the tube is infinitely long. 

In figure 10 the variation of concentration with 
distance is illustrated for a finite cathode volume 
corresponding to 7=1. The tube is assumed 
infinitely long and a^>0 for this figure. 

In figure 11 there is illustrated the variation of 
concentration along the tube for the constants 
used in eq 9f. The computation can be made by 
means of formula 9. The term n°e* outside the 
sum is referred to as the " equilibrium term". 
The way in which the successive terms in the 
series converge is illustrated by the three curves. 
The first curve is marked " equilibrium +(r )". It 
represents the effect of the term n°e^ and of the 
first term in the summation in eq 9. The second 
curve is marked " equilibrium + 0"o) + (pi)"- It- 
includes, in addition to the two previous terms, 
also the second term in the summation which is 
due to the root p x . This curve is already a fair 
approximation to the answer. The third curve 
includes the effect of the root p 2 and is marked 
a equilibrium +(r ) + (pi) + (p 2 )". It is practically 
the answer. The special value, t=0.16, was used 
for all curves in this figure. 

In figure 12 there is illustrated the superposition 
of terms for X= — 1, 7=— 8.1106, t=0.16 as in 
eq 9g. 

In this case the isotope is being pumped out of 
the cathode. The convention regarding labeling 
curves is the same as in figure 10. The convergence 
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Figure 11. 

X = 2, 7=2.6966, r=0.16. 1. Equilibrium+(r )+(p,)+(/> 2 ); 2. equilibrium 
rr ( )+(Pi)» 3. eqiiilibrtaHi+'r ^. 
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ions. For K + , the experimental value of /z is 
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Figure 12. 

X = -1,7 = -8. 1106, t=0.10. 1. Equilibrium+(r ); 2. (p^XlO 4 contribution. 

is very good. The curve lor a equilibrium-}- (r )" 
is practically tin 4 answer. The contribution due 
to the term in p x multiplied by 10 4 is shown also 
in the lower part of the figure. 

In figure 13 small X and large 7 are employed, 
X=l/5 and 7=24.022. The convergence is so 
rapid that even at r=0 only the real root need he 
considered. For 7= 1.6, ten times the r used in 
the two preceding figures, the distribution along 
the tube is not far from linear. 



'6.6X10" 4 



cm' 



The effect of convection on 



v-sec 

D h , D t due to thermal effects is not known. A 
lower limit for D h D h can be obtained from the 
relations 

which follow from statistical mechanics. Here e 
is the ionic charge and not the electronic charge. 
For a singly charged ion, e is minus the charge on 
the electron. At room temperatures for singly 
charged ions the above formulas give 

D = .026M(z>in <,m2 ; M in- C1 ^) 

\ sec v-sec/ 

For a countercurrent adjusted to balance the 
slower ion one has, assuming 



6.6X10" 4 



=3.3X10 -6 



Mz Hh~- 20Q 

a=3.3X10- 6 2? 
where E is in v/cm. Hence 
3.3X10 



enr 
V-SeC 



£= 



T.7X10" 



Ex=0.2Ex. 
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Figure 13. 

X=2, 7=24.022. Equilibrium terra+(7" ) term in both curves. 1. r=1.6; 
2. T =0. 

IV. Discussion 

1. Orders of Magnitude. 

It appears possible that the difference in 

mobilities is of the order of {^M h jMi—l)^x where /x 

is the mobility and M h , Mi are the masses of the 
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A change in £ by 1 corresponds in this case to a 
drop of 5 v along the tube. More generally, for 
the same countercurrent adjustment, and singly 
charged ions 

£1= wp=39 — times voltage drop 

Here A/x is the absolute value of the difference in 
mobilities. Thus if A m /m= 1/1,000 then 26 v are 
required to change £ by unity and to give a factor, 
e, in the concentration in equilibrium. 

The quantity r can be considered as the time in 
units D/a 2 . For K + this unit is of the order 



1.7X10" 5 
(3.3) 2 X10- 12 # 



-1.6X10 6 sec/# 2 =430 hr/E 2 



where -Bis in v/cm and the countercurrent adjust- 
ment is such as to keep one of the ions stationary. 
In the example of figure 5 one may consider r^2 as 
corresponding to the time which it takes to reach 
half equilibrium, and the time constant is in this 
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case of the order 1,000 hr/E 2 . More generally 
the natural unit of time is 



D 



(An) 2 E 2 



On account of the large uncertainty in A/z for 
uranium this unit cannot be computed with cer- 
tainty. Since A/z enters to the square, a decrease 
in A/z by a factor 10 brings about an embarrassing 
increase of a factor 100 in the time unit. For 
a small cathode volume i. e. for small y, the 
"natural" unit of time is the main thing that 
matters so that a sufficiently large A/z is then 
especially important for a small time constant. 

On the other hand, if y is large the value of r 
corresponding to half equilibrium is much greater 
than unity as is seen in the examples following 
eq 9e. According to eq 8c for a fixed X, the 
difference of mobility enters only once in the 
denominator in this case. For K + assuming 
A/x=/x/ 2 oo and taking F=20 cm 3 , A= 1 cm 2 the time 
constant is 



6X10 6 sec. 
E 



(e*-l): 



2,000 hi; 
E 



(*-l) 



It is interesting to note that if the equilibrium 
enrichment factor e x is kept fixed and if y is large 
enough to make formula 8b applicable then the 
time constant does not depend on the diffusion 
constant, D, but only on the difference of mobil- 
ity A/z. In order, however, that X be fixed for 
different D and E, one has to vary the length, L. 
For a tube 30 cm long with 90 v applied having a 
free cross-sectional area 1 cm 2 and a cathode 
volume y=20 cm 3 , one has in the absence of 
convection and for A/z//z=V2oo, A— 17.5, 7=12. 
An enormous equilibrium concentration of e 17 could 
be obtained in theory in such a case, but it would 
also take an impractical time to approach this 
equilibrium. It is doubtful, however, that D can 

/ _ cm 2 \ 

be as small as assumed ( 1.7 X 10" 5 — - J because of 

the convection currents. 

The values A=2, 7=2.6995, (used in example 
of 9f) in the absence of convection for A/z//z=72oo 
correspond to 10 v applied to the tube. If the tube 
is 30 cm long and has a free cross-sectional area 
of 1 cm 2 , the cathode compartment would have to 
be 40 cm 3 . Onehasa=3.3X10- 6 /3=l.lX10- 6 cm/ 
sec and V/A=4:0 cm so that 1 =3X 10 8 sec ~ 10 5 hr 



which is of the order of 5,000 days. If this tube 
were shortened to 3 cm keeping V and A constant 
and the same voltage then X remains unchanged 
and a is increased by a factor 10 so that T decreases 
to about 500 days. 

There is a large gain in the shortened time con- 
stant due to shortening the tube for the same volt- 
age. This is because of the increased current. 
The gain is of no real help except through the 
increased current and the detection of a difference 
in ionic mobilities appears to be easiest with a 
large current in a long tube. 

2. Equilibrium Concentration 

In this section the coefficients of diffusion will 
be decomposed into a part present in the stagnant 
liquid and a part due to spurious causes such as 
convection and stirring. 

D l =B° l +® } D h =D° h +&. 

It is supposed that the spurious effects give the 
same contribution H to B x and D h . The relations 



fXl ~kT^ h fXh ~kT^ h 



hold. In equilibrium 



where 



n t n h 
n h n, 



a niE+v ti h E J rv _ 

m-D]_ 



o(^S B ) 



The counter current is supposed to be adjusted so 
that 

v=-ii h E(l+0) 

If 0=0, then isotope h is kept stationary by the 
combined action of the countercurrent and of the 
migration velocity. However, 6 will not be sup- 
posed to be necessarily zero. One has* from the 
above equations, 



-^W+Dla+e)} 



Jf-i)l 



(m+Dl)(B+D°) 



*Prof. H. C. Urey first derived a special case of the above formula which 
corresponds to B=0, 0=0. The authors are indebted to Prof. Urey for a 
stimulating discussion in this connection. 
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For 51 = one has 

e -§ (1+*) *=f* '(9=0) 

Here at room temperatures 

TTT 

-yjr=39 x applied voltage X valence 

just as for the quantity f. It is seen from the 
above formula for a that a large 6 is advantageous 
for getting a large enrichment ratio. This fact 
should perhaps be emphasized because it indicates 
that it is better to have the countercurrent too 
strong rather than too weak if it is desired to test 
for the existence of a mobility effect. 
If B»Z>°, one has 



'' kT 



B 



-Dl^eZL^-^D* 



kT 



Vi 



(B»Z?°) 



For large © the logarithmic increment is seen to 
be a fraction Z) /j| of its value for H = 0, 0=0. 
In this limit the exact adjustment of the counter 
current is of no importance for the enrichment 
factor at equilibrium. 
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